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t"^ ■ Abstract 

O. 

The shape evolution and shape coexistence phenomena in neutron-rich nu- 
clei at A^ ~ 60, including Kr, Sr, Zr, and Mo isotopes, are studied in the 
covariant density functional theory (DFT) with the new parameter set PC- 
/\ • PKl. Pairing correlations are treated using the BCS approximation with a 

j^ ■ separable pairing force. Sharp rising in the charge radii of Sr and Zr isotopes 

at iV = 60 is observed and shown to be related to the rapid changing in 
nuclear shapes. The shape evolution is moderate in neighboring Kr and Mo 
isotopes. Similar as the results of previous Hartree-Fock-Bogogliubov (HFB) 
calculations with the Gogny force, triaxiality is observed in Mo isotopes and 
shown to be essential to reproduce quantitatively the corresponding charge 
radii. In addition, the coexistence of prolate and oblate shapes is found in 
both ^^Sr and ^°°Zr. The observed oblate and prolate minima are related 
to the low single-particle energy level density around the Fermi surfaces of 
neutron and proton respectively. Furthermore, the 5-dimensional (5D) col- 
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lective Hamiltonian determined by the calculations of the PC-PKl energy 
functional is solved for ^^Sr and ^^''Zr. The resultant excitation energy of 0^ 
state and EO transition strength p^{EO] 0^ — )■ OJi'') are in rather good agree- 
ment with the data. It is found that the lower barrier height separating the 
two competing minima along the 7 deformation in ^^''Zr gives rise to the 
larger p'^iEO; 0^ -^ Of) than that in ^^Sr. 

Keywords: 

Covariant density functional, shape evolution and shape coexistence, charge 

radii, neutron-rich Kr, Sr, Zr, Mo isotopes 

1. Introduction 

In recent decades, the evolution of nuclear shapes along isotopic and iso- 
tonic chains in neutron-rich nuclei at A^ ~ 60 has attracted many attentions. 
The sudden onset of quadrupole deformation in neutron-rich Sr and Zr iso- 
topes at the neutron number A^ = 60 is of particular interest. Such a rapid 
shape evolution has been deduced from the abrupt changing of lifetimes of 
2f states [l|, |2| as well as the quadrupole moments of rotational bands |3|. 
Besides, the excitation energies of 2^ states [4] , two- neutron separation ener- 
gies [5|, and mean-square charge radii [6[ exhibit a dramatic change between 
A^ = 58 and 60 in Sr and Zr isotopes. Very recently, the systematic of the 2^ 
states in Kr isotopes has been extended up to A^ = 60 [7|, at which nucleus, 
the energy of the first excited state drops down suddenly by ~ 400 keV. 
It indicates that the shape transition is also rather abrupt in Kr isotopes. 
However, the measured charge radii was shown to be increasing moderately 
with the neutron number at A^ = 60 |8|. 

Essential to the understanding of this dramatic shape evolution is the co- 
existence of different shapes in the two lowest O"*" states. Shape coexistence 
phenomena at low energy in Sr and Zr isotopes around A^ = 60 have been 
shown in many experimental measurements. In Ref. [9|, Jung et al. discov- 
ered two low-lying 0^ states in ^^Sr at 1229 and 1465 keV respectively. Later 
on, an extremely strong electric monopole transition of p^{EO) = 0.18 was 



observed between the first two 0"*" states 10|, lll|. The analysis of B{E2) 



and p'^{EO) values for both ^^Sr and ^°°Zr by Mach et al. indicates that these 



two nuclei have very similar structures [12[. Schussler et al. discovered a 
very low-lying 0"^ state at 215.5 keV in ^^Sr. The transition probabilities, 
the reduced EO matrix element and the observed level structure suggest the 



coexistence of a quadrupole deformed ground state and a spherical excited 



0+ state in ^^Sr [l3|. To clarify such a picture of shape coexistence, an ex- 
periment to measure the spectroscopic quadrupole moment of the 2^ state 
has been proposed [ij • 

On the theoretical side, the shape evolution around A^ = 60 has been stud- 



ied extensively with various theoretical models, including the phenomenolog- 



ical models 15, 16 




20, 21, 22 



the interacting boson model |23| . 

and the self-consistent mean-field models with 

the Gogny force 28|, |29|, |30| as well as the 



the modern shell model 

the Skyrme force 

effective relativistic Lagrangian SJJ. Most of these models have shown the 

increasing of deformations up to A^ = 60 and indeed found the competing 

prolate and oblate minima. However, the subtle balance between these two 

minima depends on the details of calculations. 

In recent years, nuclear covariant DFT has achieved great success in the 
description of ground state properties of both spherical and deformed nuclei 
all over the nuclear chart 32, l33|, |3J, |35|] . In particular, the covariant DFT 
theory with a point-coupling interaction has recently attracted more and 
more attention 36|. It shows great advantages in the extension for nuclear 
low-lying excited states by using projection technique s |37l| . generator coordi- 
nate methods 38|, l39|, |40| and collective Hamiltonian |41|]. In this framework, 
there are seve ral p opular parameter sets, including PC-Fl 42], DD-PCl 43 
and PC-PKl M 



Among these parameter sets, the PC-PKl was proposed 
very recently by fitting to observables of 60 selected spherical nuclei, including 
the binding energies, charge radii, and empirical pairing gaps. The success 
of PC-PKl has been illustrated in the description of infinite nuclear matter 
and finite nuclei for both ground-state and low-lying excited states. Further- 
more, the PC-PKl provides a good description for the isospin dependence of 
binding energy along either isotopic or isotonic chain. 

Recently, a separable pairing force with two universal parameters has 
been introduced, which was adjusted to reproduce the pairing properties of 



the Gogny force DIS in nuclear matter [45|]. The separable pairing force 



has been shown to be successful in the description of nuclear matter [45 



spherical and deformed nuclei [46|, |47|, |48|, |49|. Therefore, in this work, we 



would like to use the PC-PKl parameter set together with the separable force 
to perform a systematic calculation for the neutron-rich Kr, Sr, Zr, and Mo 
isotopes. The shape evolution and shape coexistence in this region will be 
examined. 

The theoretical framework for the relativistic point-coupling model with 



the separable pairing force is described in Sec. [21 The shape evolution in 
neutron-rich Kr, Sr, Zr, and Mo isotopes at A^ ~ 60 and shape coexistence 
phenomena in ^^Sr and ^''^Zr will be discussed in Sec. |3l Finally, a summary 
is given in Sec. HI 

2. The model 

In the covariant DFT with point-coupling interaction, the energy func- 



tional has the following form [42|, |44 



Ermf = ^ dr vl'il;k{r){-n ■ V + m)^fc(r) 
k -^ 

+ -^JTviJTv)i. + — i^y^O'Tv)/. 

+ ^-F.^F^'' - F'^doA, + e^-^j,A>^^ , (1) 

where e is the charge unit for protons and it vanishes for neutrons. The 
energy functional ([1]) contains 9 coupling constants as, ay, cxtv: Ps^ Is, 7v, 
6s, Sy and 6tv- The subscripts indicate the symmetry of the couplings: S 
stands for scalar, V for vector, and T for isovector, while the symbol refer 
to the additional distinctions: a refers to four-fermion term, 6 to derivative 
couplings, and /3 and 7 to the third- and fourth-order terms, respectively. 

The local densities and currents in the energy functional ([1]) are deter- 
mined by, 

(2) 
(3) 
(4) 



Ps{r) 


= ^vlMr)Mr), 

k 


fir) 


= )_^vlMr)rMr), 

k 


J^y(r) 


= \^4Mr)rfMr)- 



Minimizing the energy functional ([T]) with respect to ijjk, one obtains the 
Dirac equation for the single nucleons 

[j,{zd^-Vn-{m + S)]^k = 0. (5) 

The single-particle effective Hamiltonian contains local scalar S{r) and vector 
V^{r) potentials, 

5(r) = S5, K^(r) = S'^ + f.S^^, (6) 

where the nucleon isoscalar-scalar E^, isoscalar-vector S^ and isovector- 
vector S^y self-energies are given in terms of the various densities and cur- 
rents, 

S5 = asps + PspI + IspI + ^s^Ps, (7) 

S^ = avji:. + 7vU^f + SvAj>^ + eA^, (8) 

For a system with time reversal invariance, the space-like components of 
the currents and the vector potential vanish. Furthermore, one can assume 
that the nucleon single-particle states do not mix isospin, that is, the single- 
particle states are eigenstates of T3. Therefore only the third component of 
isovector potentials E^y survives. The Coulomb field Aq is determined by 
Poisson's equation. 

Pairing correlations between nucleons are treated using the BCS approxi- 
mation with a pairing force separable in momentum space, i.e., {k\V ^°\k') = 



—Gp{k)p{k') , introduced in Ref. 45| with a Gaussian ansatz p{k) = e~" ^ . 
The two parameters G and a have been adjusted to reproduce the pairing 
properties of the Gogny force DIS in nuclear matter. The obtained values 
for the parameters are G = —728 MeV ■ fm^ and a = 0.644 fm. 

In the coordinate space, the separable pairing force takes the following 
form, 

y(ri, r2, r;, r^) = GS{R- R') P(r)P(r')^ (1 - P^ , (10) 

where R = | (ri + r2) and r = ri — r2 are the center-of-mass and the relative 
coordinates respectively. -P(r) is the Fourier transform oi p{k), 

P{r) = ^e-"-'/^"'. (11) 



The pairing force has finite range, and it can preserve translational invari- 
ance due to the presence of the factor 5 (R — R'). Even though 5 (R — R') 
imphes that this force is not completely separable in coordinate space, the 
corresponding antisymmetrized pp matrix elements can be represented as a 
sum of a finite number of separable terms in the basis of a three-dimensional 
(3D) harmonic oscillator (HO): 

AfO A^^ N9 
N^=0 Ny=0 N^=0 

where N^, Ny, and A^^ are the quantum numbers of the corresponding one- 
dimensional (ID) HO in the center-of-mass frame. The summations over 
Nx, Ny, and A^^ are restricted to finite terms with cutoffs N^,Ny, and N^ 
respectively. The convergence with respect of the cutoffs has to be checked 
in calculations. V^^ " "^ represents the single-particle matrix element in the 
3D HO basis. In this case, the pairing field can be written as a sum of a 
finite number of separable terms 

n2 no atO 

^"^ = G E E E (yT-"rp....... ds) 

N^=0 Ny=0 N:,=0 



with the coefficients 



PN.Ny.. = E Vr^-^n,,, (14) 

75>0 



where k^j is the matrix element of pairing tensor. The expression of V 



has been derived in Ref. 48 



76 



In the BCS approximation, the pairing gap A^ for each single-particle 
state -0^ is finally determined as follows, 

Ak = ^A^-pFkaFk(3, (15) 

where F^a is the expansion coefficient for the large component in Dirac spinor 
i/jk on the 3D HO basis. The resultant pairing energy is given by 

n2 no no 

Epi,n = *^ 2^ 2^ 2^ {PN^NyN.yPN^NyN,- (16) 

N^=0 Ny=0 N^=0 

6 



The center-of-inass correction to the energy is considered microscopically 
with both the direct and exchange terms, 

^„.., = -^. (17) 

where m is the mass of nucleons. A is mass number and Pcm = Xli Pi i^ ^^^ 
total momentum in the cm. frame. 

The total nuclear energy is determined by 

-E'tot = -E'rMF + -E-pair + -E'c.m. (18) 

The potential energy surface (PES) in the plane of deformation variables is 
obtained by imposing a quadratic constraint on the mass quadrupole mo- 
ments 

(i/)+5^C2,((g2M)-92,)' (19) 

At=0,2 

where {H) is the total energy, and {Q2^l) denotes the expectation value of 
the mass quadrupole operator: 

Q20 = 2z'-x^-y^ (20) 

Q22 = x^-v" (21) 

Here q^^ is the constrained value of the quadrupole moments, and 6*2,^ the 



corresponding stiffness constant (50 



3. Results and discussion 



The new parametrization PC-PKl [4J] and the separable pairing force [45 



are adopted in the particle-hole channel and the particle-particle channel 
respectively. Parity, D^ symmetry, and time-reversal invariance are imposed. 
The Dirac equation is solved by expanding in the basis of eigenfunctions of 
a 3DH0 in Cartesian coordinate with 12 major shells, which are found to be 
sufficient to obtain a reasonably converged mean- field PES. 

3.1. Shape evolution in neutron-rich Sr isotopes 

Figure [U displays the PESs of even-even 88-i04gj, -^^ p_^ plane, normalized 
to the total energy of absolute minimum. The energy difference between 
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Figure 1: (Color online) The potential energy surfaces of the even-even 88-i06gj. jgotopes 
in the /3-7 plane from the constrained relativistic mean-field (RMF) plus BCS calculations. 
All energies are normalized with respect to the tot energy of the absolute minimum. The 
energy difference between neighboring contour lines is 0.5 McV. 



4.6 



4.5 



4.4 



4.3 



4.2 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


—A— sph. ,# [ 


-•--pro. ,•' 


- -T-obl. T ,•' 


•' A -X'^- 


J ^--y ■ 


7/* .^' 


/r > 




Sr 


" 


M ' J 

1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 . 1 



50 52 54 56 58 60 62 64 66 68 

N 



Figure 2: (Color online) The evolution of the nuclear charge radii in Sr isotopes. The 
calculated values corresponding to the spherical (up triangles) , prolate (circles) , and oblate 
(down triangles) local minima in the PESs (Fig. [T]) are plotted as functions of neutron 
number. The squares with error bars denote the experimental data [51| . 



neighboring contour lines is 0.5 MeV. The PESs in Fig. [T]show a clear pic- 
ture for the evolution of shapes in ^^"^''^Sr. Starting from a well spherical 
shape of ^^Sr, the spherical (global) minima in ^'''^^Sr become soft against 
the distortion towards oblate shape. In the mean-time, the prolate (second) 
minimum comes down. When the neutron number increases from A^ = 56 
to A^ = 60, the global minimum is shifted to the oblate side with large de- 
formation. Meanwhile, the prolate minimum becomes deep and competing 
with the oblate minimum at A^ = 60. A triaxial barrier with the height 
~ 2.22 MeV separates these two competing minima in ^^Sr. Beyond the 
A^ = 60, the structure of the energy maps is stable, that is, a soft oblate 
minimum against the distortion towards spherical shape coexists with a well 
prolate one. 

The evolution of the nuclear charge radii in neutron-rich Sr isotopes can 
be seen in Fig. [21 where the calculated charge radii corresponding to the 
spherical, prolate, and oblate local minima in the PESs of even-even ss-io^g^, 
(c.f. Fig. [1]) are plotted as functions of neutron number. For A^ > 54, the 
charge radii of spherical shapes are given as well for comparison. It is shown 
that the charge radii of the spherical and prolate shapes increase smoothly in 
the similar slop with the neutron number. The difference in the two charge 
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radii is about 0.12 fm, originating from the effect of prolate deformation. 
Moreover, it means that the deformation of prolate minimum is nearly the 
same when the neutron number increases from A^ = 60 to A^ = 68. On the 
contrary, the charge radius corresponding to the oblate minimum changes 
rapidly with the neutron number. In particular, a sudden rising of charge 
radius from A^ = 58 to A^ = 60 and a sudden dropping from A^ = 60 to 
A^ = 62 are due to the increasing and decreasing of the oblate deformation 
from j3 = —0.25 to /? = —0.35 and back to /3 = —0.2. Moreover, it is shown 
in Fig. [2] that the charge radii of prolate and oblate minima in ^^Sr are 
similar due to the similar size of quadrupole deformation. Comparing with 
the available data for charge radii, one can draw a shape evolution picture for 
the ground states of even-even *^~^°°Sr, namely, from spherical shape (^^Sr) 
to more oblate shape (^'^Sr), oblate and prolate coexistence (^*Sr) and finally 
more prolate shape (^"'^Sr). 

3.2. Shape evolution in neutron-rich Kr, Zr, and Mo isotopes 

The PESs of neutron-rich Kr, Zr, and Mo isotopes are shown in Figs. |3]|Sl 
respectively. In comparison with the shape evolution picture of Sr isotopes, 
the main difference is found in the evolution of prolate minimum in Kr iso- 
topes, where the PESs are much softer and the prolate minima are not well 
developed. 

The shape evolution picture of Zr isotopes is very similar as that in Sr 
isotopes, except the barrier height separating the prolate and oblate minima. 
In contrary with the case in Sr isotopes, the prolate and oblate minima are 
always connected through triaxial distortion with near-zero barrier height, 
in particular, for ^°°Zr with shape coexistence phenomenon. 

For Mo isotopes, the shape evolution picture is similar as that in Zr 
isotopes. The evident difference is the occurrence of triaxial minima in the 
Mo isotopes with neutron number from A^ = 58 to A^ = 68. 

Very recently, a global study of nuclear low-lying states based on the non- 
relativistic Hartree-Fock-Bogoliubov framework with the Gogny force have 
been done. The corresponding potential energy surfaces and other observ- 



ables are given in Ref. [29|. Based on the same framework, Rodriguez- 
Guzman et al. have examined in detail the shapes evolution of nuclear 
ground-state in neutron-rich Sr, Zr, and Mo isotopes, including both even- 



even and odd- A nuclei [30|. The trend of shape evolution is similar as our 



result based on the covariant density functional. However, the transition at 
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Figure 3: (Color online) Same as the Fig.[TJ but for the isotopes ^^ -'^^'^Kr. 
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Figure 4: (Color online) Same as the Fig. [1] but for the isotopes ^° ^"^Zr 
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Figure 5: (Color online) Same as the Fig.[Tl but for the isotopes ^^ '^^"Mo. 
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Figure 6: (Color online) Same as the Figl^J but for Kr isotopes. 

A^ = 60 in our calculations is a little slower along the isotopic chain, and 
more rapid along the isotonic chain for the A^ = 60 isotones. 

The charge radii of Kr, Zr, and Mo isotopes are plotted in Figs. |6]|8l 
respectively. The sharp transition is also observed in Zr isotopes, which 
also indicates the rapid change in the nuclear shapes. On the contrary, 
the charge radii in Kr and Mo isotopes increase smoothly with the neutron 
number. Similar as the calculation results of the Gogny force in Ref. 30|, 
the triaxiality is shown to be essential to reproduce qualitatively the charge 
radii in Mo isotopes. 



3.3. Covariant density functional based 5D collective Hamiltonian analysis 
of shape coexistence in ^®Sr and -"^"^Zr 

The coexistence of prolate and oblate shapes observed in ^^Sr and ^''^Zr 
will be studied in more detail with the 5D collective Hamiltonian determined 
by the constrained self-consistent RMF plus BCS calculations. The details 
about the covariant density functional based 5D collective Hamiltonian can 
be found in Ref. 



48 



In Fig. [9] we plot the total energies as functions of axial deformation /3 for 
^*Sr and ^°°Zr. The inset displays the PEC corresponding to the projections 
on the 7 deformation, that is, the minimum for each 7 deformation on the 
PES in the /3-7 plane (c.f. Figs. [UandHj). In both nuclei, the coexisting 
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Figure 7: (Color online) Same as the Figl^l but for Zr isotopes. 
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Figure 8: (Color online) Same as the Figl^l but for Mo isotopes. 
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Figure 9: (Color online) The total energies of ^*Sr and ^™Zr as functions of axial defor- 
mation /?. All energies are normalized with respect to the total energy of the absolute 
minimum. The inset displays the PEC corresponding to the projections on the 7 defor- 
mation, that is, the minimum for each 7 on the PES in the /3-7 plane. 



prolate and oblate minima with very closed binding energies are observed, 
which are separated by certain barriers. In ^^Sr, the spherical barrier height 
is ~ 3.4 MeV. After considering the 7 degree of freedom, this barrier height 
is lowered down to 2.2 MeV. In ^'^''Zr, the barrier height is much smaller 
with the size ~ 0.5 MeV if the 7 deformation is considered. 



Table 1: The calculated excitation energies (in MeV) of Oj states and EO transition 
stren gth s p^jEO; Oj -^ 0^) x 10^ in ^^Sr and ^""Zr, in comparison with the corresponding 
data ^, [5I 





98Sr 


ioo2r 




Cal. Exp. 


Cal. Exp. 


E(0+)(MeV) 
p2(E0;0+^0+) X 103 


0.216 0.215 
116.841 51(5) 


0.468 0.331 
150.321 108(19) 



The excitation energy of the 0^ state and the EO transition strength 
p^(£'0; 0^ -^ 0^) between the 0^ and Of states are two key quantities in the 
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study of shape coexistence, 

(0^IE.efcr2|0+) 



p2(E0;0+->0+ 



eRl 



(22) 



where Rq ~ 1.2^4^/^ fm. The p^{EO; 0^ — )■ Of) is related to the change in the 
root mean-square charge radius of the nucleus between the Of and 0^ states, 
and therefore carries important information about the change in deformation 
and the overlap of the wave functions. 

In Tab. [H we list the calculated excitation energies of 0^ states and EO 
transition strengths p^{EO; Of — )■ 0^) in ^^Sr and ^°°Zr from the solution of 
5D collective Hamiltonian based on the energy functional PC-PKl plus the 



separable pairing force. The experimental data [52|, |53| are also shown for 
comparison. The existence of very low-lying 0^ state is often used as a strong 
signal for the shape coexistence. As expected, the calculated excitation en- 
ergies of 0^ states in both nuclei are predicted in very low values, that is, 
0.216 MeV for ^^Sr and 0.468 MeV for ^°°Zr, which are also very close to the 
data. Although the experimental EO transition strengths p^(i?0;0^ — )■ 0^) 
are overestimated by the collective Hamiltonian based on PC-PKl functional, 
they are all typically large, again confirming the shape coexistence phenom- 
ena in these two A^ = 60 isotones. 

The mixing between the 0^ and 0^ states can be further understood from 
the distribution of the wave functions of the 0^ and 0^ states. Figure dU] 
displays the probability density distribution of Of and 0^ states in /3-7 plane 
for ^^Sr and ^"'^Zr. Due to the hight triaxial barrier (c.f. Fig. [1]), two peaks 
corresponding to the coexisting prolate and oblate shapes are observed in 
both O]*" and 0^ states in ^^Sr. However, in ^°°Zr, the probability density 
of the O^i*' state is almost uniformly distributed along the 7 deformation, 
connecting the prolate and oblate shapes. Since there is one node in the 
probability distribution of 0^ state, the resultant p^i^EO-jOf — ;■ 0^) is the 
consequence of concelation from the probability distributions of prolate and 
oblate parts. This cancelation is larger in ^^Sr than that in ^'^'^Zr. As a result, 
the obtained p^{EO\ 0^ — > 0^) in ^^Sr is much smaller than the value in ^''^Zr, 
as shown in Tab. [T] 

The observed shape coexistence phenomenon can be understood from the 
distribution of single-nucleon levels. In Fig. [TTl we plot the neutron and pro- 
ton single-particle energy levels in ^^Sr as functions of the axial deformation 
parameter j3. The thick dotted curves denote the position of the correspond- 
ing Fermi levels. It is shown that the neutron Fermi level goes across the 
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Figure 10: (Color online) Probability density distribution in the (3-^ plane for the Q^ and 



Oj states of ^^Sr and ^^Zr. 
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Figure 11: (Color online) Neutron and proton single-particle levels for ^*Sr as functions 
of the axial deformation parameter /?. The thick dotted curves denote the corresponding 
Fermi levels. 



deformation region of low level density with —0.4 < /3 < —0.15, giving rise 
to the oblate minimum. On the other hand, the proton Fermi level locates in 
the middle of the energy gap around (3 ~ 0.45, which gives rise to the prolate 
minimum in ^^Sr. 

4. Summary 

In summary, the triaxial relativistic mean-field plus BCS model with a 
point-coupling interaction in the particle-hole channel and a separable pairing 
force in the particle-particle channel has been established and applied to 
study the shape evolution and shape coexistence phenomena in neutron-rich 
A ~ 100 nuclei, including Kr, Sr, Zr, and Mo isotopes using the newly 
parameterized PC-PKl energy functional. The evolution of potential energy 
surfaces and charge radii with the neutron number in each isotopes have 
been presented. Sharp rising in the charge radii of Sr and Zr isotopes at 
A^ = 60 has been observed and shown to be related to the rapid changing in 
the nuclear shape. This dramatic evolution of charge radii is smoothed out 
in Mo isotopes due to the occurrence of triaxial minimum, which is similar as 
the results of Hartree-Fock-Bogogliubov calculations with the Gogny force. 
In particular, the triaxiality has been shown to be essential to reproduce 
quantitatively the charge radii of Mo isotopes. 

The coexistence of prolate and oblate shapes has been observed in ^^Sr 
and ^^^Zt. However, the barrier height separating the coexisting minima 
along the 7 deformation in ^°°Zr has been shown much lower than that in ^^Sr. 
The observed oblate minimum and prolate minimum are related to the low 
single-particle energy level density around the Fermi surfaces of neutron and 
proton respectively. Furthermore, the 5D collective Hamiltonian determined 
by the calculations of the PC-PKl energy functional has been constructed 
and solved for ^^Sr and ^°°Zr. The resultant excitation energy of 03" state 
and EO transition strength p^(i?0;0^ — ?■ 0^) are in rather good agreement 
with the data. It has been found that the lower barrier height in ^°°Zr gives 
rise to the larger p^i^EO; 0^ — t- 0^) than that in ^^Sr. 
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